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Abstract. Using Griffiths and Lieb-Simon type inequalities, it is shown that 
the two-point function of ferromagnetic spin models with N components in 
one dimension decays like the interaction J(n) ~ n~ 7 provided that 1 < N < 
4 and T > T c . 



1 Introduction and Main Result 



As is well known, classical spin models with a continuous symmetry in two 
dimensions lead to scale invariant field theories with the nonlinear sigma- 
model action 

J rf 2 *(Vs(x)) 2 , (1.1) 

where s is a unit spin with iV components. The resulting behaviour distin- 
guishes between an abelian, N = 2 plane rotor, and a non-abelian symme- 
try group, N > 3. In the latter case the two-point correlation, g(x) =< 
s(x) • s(0) > decays exponentially at any finite temperature with a finite 
correlation length £(T) ~ exp[27r/(iV — 2)T] |T], QL On the other hand for 
an XY-symmetry, the exponential decay holds only above the Kosterlitz- 
Thouless critical temperature Tkt- The phase at low temperatures exhibits 
power-law decay g(x) ~ |x| _T '^ with a continuously varying exponent rj(T) 
[||]. Qualitatively this behaviour can be understood through the spin- wave 
approximation 

(Vs) 2 « (V^) 2 (1.2) 

with s = (cos (p, sin (p), neglecting the periodicity of the phase variable (p. 
Since the approximate action is Gaussian, the correlation function <?(x) = 
9ft < exp[i(ip(x) — (p(0))] > can be calculated easily, yielding the power law 
decay |x| _,? ^ with rj = T/2n. The vortex excitations lead, at low T, only to 
a finite renormalisation of r\ P|. 

In one-dimensional models the spatial dimension can be mimicked by a 
long range interaction with a decay as |n| -7 , n being a point on the one- 
dimensional lattice. One notices that for 7 = 2 the action is again scale 
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invariant. For this marginal case, it was conjectured early by Thouless M 
that the N = 1 Ising model has a spontaneous magnetization m* below 
a nonzero critical temperature T c . The spontaneous magnetization jumps 
to a finite value at T c , yet the transition is continuous. This conclusion 
was confirmed through an analysis of the equivalent Kondo problem ||, ^| 
and later proven rigorously 0, §]. A renormalization group calculation for 
long range spin models in one dimension both for iV = 1 and a continuous 
symmetry was performed by Kosterlitz || . Within a one loop calculation he 
showed that there is always a low temperature spontaneous magnetization 
provided 1<7<2(7>1 being necessary to have an extensive free energy). 
In the marginal case 7 = 2 and for N > 2, the associated beta function 
vanishes quadratically near the trivial fixed point T = 0. This indicates that 
T c = for N > 2 and an exponential behaviour x{T) ~ exp[2ir 2 /(N — 1)T] 
for T — > of the susceptibility. Due to the power law interaction there can 
be no finite correlation length, however. 

In our present note we discuss the long distance behaviour of the two- 
point function g(n). In the phase where m* = it is shown that if g(n) has 
a power law decay at all, it is necessarily equal to that of the interaction. 
In particular, a spin wave approximation is qualitatively incorrect for long 
range models even at very low temperatures. Moreover, in the XY-case 
with marginal 7 = 2, our rigorous bounds rule out the appearance of a low 
temperature phase with a continuously varying exponent f]{T) and infinite 
susceptibility, which has been claimed in the literature on the basis of spin 



wave theory and Monte Carlo simulations 1 10, 11, 12]. Our result is of direct 



relevance to the problem of strong tunneling in the so called single electron 



box, showing that the Coulomb blockade at zero temperature is not destroyed 



even for large conductance |L3 



To be more precise, we consider the spin Hamiltonian 

H = -ij^ J(m-n)s m -s n , (1.3) 

m,n 

where the couplings are ferromagnetic, J(n) > 0, and decay as 

J(n) = \n\-^ (1.4) 

for |n| — > oo. As before the two-point function is defined by g(n) =< s h -Sq > 
in the infinite volume limit with free boundary conditions. Then 

lim g{n) = m* 2 , (1-5) 

where m* is the spontaneous magnetization with the standard convention 
that m* = for T > T c and m* > for T < T c . We define the scaling 
exponent rj by 

g(n) -m* 2 = |n|~" (1.6) 
for large n. The magnetic susceptibility is given by 

X=^J2(9(n)~m* 2 ). (1.7) 



If r) < 1, then x = oo. 

The qualitative phase diagram for such feromagnetic models is rather well 
understood. For 7 > 2 one has m* = and for 7 < 2 m* > at sufficiently 
low temperatures. In the marginal case, 7 = 2, the number of components 
becomes relevant. Whereas for the Ising model, N — 1, m* > at low T [pj], 



for N > 2 Simon proves that m* = at any finite T. The decay of the 
two-point function has been studied on a rigorous level mostly for the Ising 
model with particular attention to the marginal case 7 = 2 fl4]| : For T > T c 
one has rj = 2. At T c m* jumps to a non-zero value, the Thouless effect, and 
77 jumps to zero. Below T c , rj(T) increases with decreasing T and locks to its 
high temperature value 77 = 2 at some critical value T* < T c . Although rj 
varies continuously, the overall behavior is obviously quite distinct from the 
standard Kosterlitz-Thouless scenario in the short range d = 2, N = 2 case. 

Here we show that if N — 1,2,3,4, if m* = 0, and if g(n) is known to 
have some decay already, then r] = 7. A lower bound of this form is known 
from Griffiths second inequality for arbitrary iV |1| . A corresponding upper 



bound is slightly more involved. It uses a Lieb-Simon type inequality |fL7| , [TS 
which relies on Gaussian domination of the four-point function. Although 
this is expected to hold in general, it has been proved only for iV = 1, 2, 3, 4 
components |19| 



In the following section we give the details of the argument. In fact, it 
would be of interest to have a numerical solution of the nonlinear integral 
equation (2.8), which could be used as a sharp test of Monte-Carlo simula- 
tions. 



2 Bounds on the two-point function 

We consider ferromagnetic spin models with iV components in one space 
dimension with Hamiltonian (1.3). Here s ra is the iV component spin at 
lattice site n, n integer, with |s n | = 1. The couplings J satisfy J(n) > 
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0, J(n) = J{—n) and have the asymptotic decay (1.4). To have an extensive 
free energy we require 7 > 1. The equilibrium distribution in finite volume 
[—£, ...,£} is given by 

e 

Z- 1 cxp[-(3H] H 5(\s n \ - l)d N s n . (2.1) 
n=-e 

We choose free boundary conditions, i.e. s n = for n outside [—£, ...,£] and 
denote the corresponding expectation by < ■ >g. The two-point function in 
the infinite volume limit I — > 00 is then defined by 

g(m - n) =< s m ■ s n >= lim < s m • s n > £ > 0. (2.2) 

If m* — 0, g is independent of the boundary conditions. 

To discuss the asymptotic decay of g we first note that by Griffiths second 
inequality g is increasing in the couplings. Thus 

g(n) > J <J(|s | - l)d N s 5(\s n \ - l)d N s n exp[pJ(n)s ■ s n ]s • s n (2.3) 

which proves that g(n) cannot decrease faster than the couplings J(n). 

The upper bound for g is slightly more complicated and uses the well 
known Lieb-Simon type inequality. We define 

Al = {u,v\ either \u\ < L, \v\ > L or \u\ > L, \v\ < L} (2.4) 

and split the Hamiltonian as 

H\ — Hi + XH 2 , 

£ 



H-l = ~ Jim- n)s m -s n , 

m,n=—£ m,n&h c L 
1 1 

-- Yl J(m-n)s m -s n . (2.5) 



2 

m,n=—£ m,n6Ai 



Differentiating with respect to A we obtain 

f l d 

< s m ' s n >£—< S m • S n >£A=0 + / d\ — < S m ■ S n >£ \ 

Jo rf A 

dX - ^2 J(u- v) 

u,v=—£ u,v€A L 

(< ( s m - s n)( s « ' s v) >£,X — < s m ■ s n >£,X< s u ' s v >l,x)- (2-6) 

We choose \m\ < L, \n\ > L. Then the first term in (2.6) vanishes. For the 
second term we use the Gaussian domination valid for N = 1, 2, 3, 4 



< ( s m ' s n)( s « ' s v) >^,A<< s m ' s n >£,X< s u ' s v > £,X 
+~jy{< S m • S u >£ ; a< S m ■ > e x + < s m " s v >£,X< s n ' s u >l,x) (2.7) 

and set A = 1 because < s m ■ s n >t t \ is increasing in A. Finally we take 
£ — ► oo and arrive at 

</(n) < (/3/iV) E ^ J ( w - " n ) ( 2 ' 8 ) 

|u|<L |y|>L 

for \n\ > L. 

The integral inequality (2.8) is studied in [p0[| . In essence, one splits the 
v-sum into terms with \v\ < |ti/2| and those with |i>| > |n/2|. This yields for 
\n\ > L , L fixed, 

g{n) < c'\n\-^ + a{n)g(n/2) (2.9) 
with a(n) — > as \n\ — > oo. Iterating (2.9) results in a bound as 

#(n)<c|n|~ 7 . (2.10) 



The precise conditions, cf. p0 |, Lemma 5.4, for the validity of (2.10) are (i) 



for 7 > 2 it is required that lim n _ +00 g(ri) = 0, which we know already from 
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m* = 0, (ii) for 7 = 2 it is required that g(n) < ci(l + log(l + M)) -1 with 
a suitable constant c\ depending on the prefactor of J, (iii) for 1 < 7 < 2 it 
is required that g{n) < c 2 (l + M) 7-2 with a suitable constant c 2 depending 
on the prefactor of J. We conclude that under the stated conditions on g(n) 
and if N = 1, 2, 3, 4, m* = 0, then 



g{n) ~ const. |n| 7 (2-H) 

for large \n\. 

For Ising spins the bounds (2.3), (2.10) have recently been sharpened [pi 



such as to determine also the prefactor in (2.11). Generalizing to the present 
case, we conjecture that 

lim -— - < s • s n >= ^N X 2 - (2.12) 



The proof in [^TJ uses the FK and percolation representation for the lower 
bound and the random current representation for the upper bound, which 
unfortunately are special to N — 1. 
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